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Conical  convergent-divergent  (CCD)  nozzles  represent  an  important  category  of  supersonic 
jet-engine  nozzles  which  require  variable  throat  areas  and  variable  exit  areas  to  adapt  to  a  range  of 
operating  conditions.  CCD  nozzles  with  design  Mach  numbers  of  1.3,  1.5,  and  1.65  are  examined 
experimentally  over  a  range  of  fully  expanded  Mach  numbers  from  1.22  to  1.71.  The  characteristics 
of  the  flow  and  acoustic  fields  from  these  nozzles  are  explored.  Shadowgraph,  Particle  Image 
Velocimetry,  far-field  and  near-field  acoustic  surveys  are  presented.  Results  of  a  Monotonically 
Integrated  Large  Eddy  Simulation  are  presented  for  the  Mach  1.5  nozzle  at  an  underexpanded 
condition.  The  agreement  between  simulations  and  measurements  is  excellent.  It  is  shown  that 
these  nozzles  differ  from  traditional  smoothly  contoured  method-of-characteristics  nozzles  in  that 
they  never  achieve  a  shock  free  condition.  Furthermore  it  is  shown  that  these  nozzles  produce 
a  “double  diamond”  pattern  in  which  two  sets  of  shock  diamonds  are  generated  with  an  axial 
displacement  between  them.  The  cause  of  this  phenomenon  is  explored.  It  is  further  shown  that  as 
a  consequence  they  are  never  free  from  shock-associated  noise  even  when  operated  at  perfect 
expansion.  In  spite  of  this  difference,  it  is  found  that  CCD  nozzles  behave  like  traditional 
convergent-divergent  nozzles  in  that  they  produce  the  same  shock-cell  size,  broadband 
shock-associated  noise  peak  frequency,  and  screech  frequency  as  traditional  convergent-divergent 
nozzles.  The  apparent  source  regions  for  mixing  noise,  broadband  shock  associated  noise  and 
screech  are  all  similar  to  those  from  traditional  convergent-divergent  nozzles.  ©  2011  American 
Institute  of  Physics,  [doi:  10. 1063/1. 3657824] 


I.  INTRODUCTION 

Since  the  early  twentieth  century,  it  has  been  recognized 
that  convergent-divergent  nozzles,  when  operated  off-design, 
generate  semi-periodic  shock  diamonds  as  shocks  and 
Prandtl -Meyer  fans  alternate  in  the  jet  plume.1  For  traditional 
smoothly  contoured  convergent-divergent  (C-D)  nozzles, 
designed  for  isentropic  expansion  and  parallel  exit  flow, 
matching  the  exit  pressure  to  the  ambient  pressure  produces  a 
supersonic  jet  which  is  shock  free.  Such  a  condition,  with 
pressures  matched  at  the  exit,  is  referred  to  as  the  nozzle’s 
design  condition. 

In  order  to  attain  a  shock  free  design  condition,  tradi¬ 
tional  convergent-divergent  nozzles  (C-D  nozzles)  employed 
in  most  laboratory  studies  have  special  design  features.  The 
nozzle  throats  are  smooth  with  large-radius  curves.  Isen¬ 
tropic  expansion  in  the  divergent  section  is  achieved  by 
designing  the  nozzle  shape  using  the  method  of  characteris¬ 
tics  or  by  numerical  simulation.  The  contours  of  these  tradi¬ 
tional  nozzles  have  zero  slope  at  the  exit  so  that  they 
produce  parallel  exit  flows.  Such  contours  present  a  practical 
problem  for  the  designer  of  a  high-performance  jet  engine 
exhaust  nozzle.  Since  aircraft  must  operate  over  a  range  of 
atmospheric  conditions  and  power  levels,  a  variable- 
geometry  nozzle  is  generally  required.  The  throat  area  must 
be  changed  to  attain  the  desired  mass  flows  and  the  area  ratio 
between  throat  and  exit  must  be  changed  to  adapt  to  condi¬ 
tions  and  thrust  requirements. 


While  it  is  possible  to  design  a  variable  geometry  nozzle 
with  smooth  contours,  such  a  nozzle  would  be  complex  and 
heavy.  What  is  commonly  done  in  high-performance  jet  air¬ 
craft  exhausts  is  to  employ  a  variable  geometry  nozzle  with 
flat  plates  in  the  convergent  section  and  flat  plates  in  the  di¬ 
vergent  section  with  seals  between  so  that  the  plates  make  a 
cross  section  that  is  polyhedral  and  approximates  a  circle. 
The  plates  are  hinged  between  the  convergent  and  divergent 
sections  making  a  relatively  sharp  corner  at  the  throat.  Such 
a  nozzle  trades  a  reduction  in  nozzle  efficiency  for  a  reduc¬ 
tion  in  weight  and  complexity.  The  conical  C-D  (CCD)  noz¬ 
zles  employed  in  this  study  match  the  stream-wise  contours 
of  the  real  engine  nozzles.  The  cross  sections  are  circular 
and  thus  the  convergent  and  divergent  sections  are  conical. 

Like  the  variable  geometry  nozzles  used  in  actual  prac¬ 
tice,  these  conical  nozzles  do  not  approach  a  cylindrical 
shape  at  the  exit,  so  unless  the  flow  is  massively  separated 
they  can  never  emit  parallel  flow.  Turning  of  the  flow  at  the 
exit  will  be  necessary  and  this  will  produce  shocks  or 
Prandtl-Meyer  waves  in  the  jet  even  at  the  design  condition 
when  there  is  no  pressure  mismatch.  It  will  be  shown  that 
these  nozzles  are  not  shock  free  at  or  near  their  design  condi¬ 
tions.  It  will  also  be  shown  that  these  nozzles  generate  a  sec¬ 
ond  set  of  shock  diamonds  superimposed  on  the  first, 
generating  a  double-diamond  pattern.  This  second  set  of 
shock  diamonds  has  its  origin  at  the  sharp  nozzle  throat.  The 
presence  or  absence  of  shocks  in  the  jet  plume  has  a  signifi¬ 
cant  effect  on  the  acoustic  emissions  of  the  jet. 
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A  convergent-divergent  nozzle,  whether  conical  or  tradi¬ 
tional,  may  be  characterized  by  a  design  Mach  number,  Md, 
which  is  the  Mach  number  at  the  nozzle  exit  corresponding 
to  the  ratio  of  areas  at  the  nozzle  throat  and  nozzle  exit  by 


M] 


Lr 


(y+i)/(v— i) 


where  A*  and  Ae  are  the  throat  and  exit  areas  respectively 
and  y  is  the  ratio  of  specific  heats.  Only  sonic  or  supersonic 
solutions  are  of  interest  so,  Md  >  1 .  The  operating  condition 
of  an  exhaust  nozzle  may  be  characterized  by  a  jet  Mach 
number,  Mj,  which  corresponds  to  the  pressure  ratio  between 
the  total  pressure  supplied  to  the  nozzle  and  the  ambient 
static  pressure  outside  the  nozzle  by 


Po 

Poo 


y/b-i) 


where  pQ  and  p.x  are  the  total  pressure  and  ambient  pressure, 
respectively,  and  y  is  again  the  ratio  of  specific  heats.  The 
non-negative  Mach  number  solution  is  used.  The  nozzle  is 
then  said  to  be  overexpanded  when  Mj  <  Md,  on-design 
when  Mj  =  Md,  and  underexpanded  when  Mj  >  Md. 


A.  Acoustics  of  shock-containing  jets 

The  noise  generated  by  shock-containing  supersonic  jets 
can  be  divided  into  three  dominant  components:  screech, 
broadband  shock-associated  noise,  and  turbulent  mixing 
noise.  The  turbulent  mixing  noise  is  not  unique  to  supersonic 
jets  and  has  been  dealt  with  extensively  in  the  literature  of 
subsonic  jets  classically2-5  and  more  recently  in  terms  of 
large-scale  structure  analysis.6-8 

Broadband  shock-associated  noise  may  be  identified  by  a 
spectral  peak  whose  frequency  is  a  function  of  the  direction  of 
radiation.  The  shock-associated  noise  peak  frequency 
increases  in  the  direction  of  the  jet  flow.  The  frequency  de¬ 
pendence  on  the  angle  of  observation  was  first  explained  as 
constructive  interference  of  emissions  from  a  phased  array  of 
correlated  sources  caused  by  large-scale  structures  passing 
through  the  regularly  spaced  shock  diamonds.6  A  mathemati¬ 
cally  more  complete  modeling  of  the  large-scale  structures  in 
the  shear  layer  reproduces  the  frequency  dependence  of  the 
earlier  theory,  but  also  accounts  for  the  broadening  of  the 
spectral  peak  with  increasing  inlet  angle  and  the  strong  reduc¬ 
tion  in  intensity  of  higher-order  peaks.10  Broadband  shock- 
associated  noise  is  radiated  most  strongly  in  the  upstream 
direction  and  grows  stronger  with  (Md2  -  M  2)2  up  to  a  maxi¬ 
mum  amplitude  at  the  point  at  which  Mach  disks  form  in  the 
jet. 1 1  The  majority  of  broadband  shock-associated  noise  has 
been  found  to  occur  in  the  middle12  or  in  the  later  part 1 1  of  the 
region  in  which  shock  cells  exist. 

When  first  encountered,  screech  was  observed  to  be 
caused  by  an  interaction  between  shear  layer  turbulence  and 
the  shock  waves  in  the  jet.14  The  near-field  acoustic  pressure 
distribution  for  screech  has  been  shown  to  have  a  lobed  char¬ 
acter  with  regions  of  high  intensity  corresponding  to  the  shock 
reflections  from  the  shear  layer.15  The  upstream -running 


waves  in  a  screeching  jet  have  been  found  to  travel  at  the 
speed  of  sound  and  so  it  may  be  inferred  that  they  are  them¬ 
selves  acoustic  waves.16  Later,  the  importance  of  instability 
and  large-scale  coherent  structures  in  turbulence  have  been 
recognized.  Instability  waves  interacting  with  the  shock  struc¬ 
ture  and  with  sound  waves  outside  the  jet  have  been  shown  to 
be  essential  to  the  generation  of  screech. 1 7  The  flow  and  the 
acoustics  of  circular  and  non-circular  supersonic  jets  have 
been  compared.18 

The  primary  peak  frequency  of  screech  is  never  greater 
than  the  broadband  shock-associated  noise  peak.  In  fact,  it 
matches  the  primary  peak  frequency  for  broadband  shock- 
associated  noise  for  an  observer  directly  upstream  (i/^  =  0°). 
This  observation  has  led  to  the  theory  that  the  broadband 
shock-associated  noise  is  itself  the  mechanism  of  the  upstream 
running  portion  of  the  feedback  loop  producing  screech. 1 9 

Until  the  last  few  years,  nearly  all  published  results 
relating  to  shock  noise  from  supersonic  jets  from  C-D  noz¬ 
zles  were  obtained  using  traditional,  smoothly  contoured 
nozzles.  The  one  exception  reported  far-field  overall  sound 
pressure  level  (OASPL)  for  a  nozzle  described  as  “an  FI 5 
nozzle”  which  showed  elevated  noise  near  the  design  condi¬ 
tion  as  compared  to  a  traditional  C-D  nozzle.20  Recent  inter¬ 
est  in  noise  from  military  aircraft  has  motivated  a  number  of 
studies  of  the  acoustics  of  jets  from  the  kinds  of  variable- 
geometry  nozzles  used  on  military  jets.  A  high-fidelity  noz¬ 
zle  with  the  facets  associated  with  the  flat  panels  and  seals  in 
a  real  engine  nozzle  has  been  tested  and  far-field  acoustics 
presented  in  one-third  octave  bands  with  and  without  several 
noise  reducing  strategies.21 

The  design  of  the  nozzle  used  by  Martens  and  Haber21 
was  used  by  the  present  authors  as  a  guide  in  designing  the 
nozzle  used  in  this  study.  Initial  narrow-band  far-field  meas¬ 
urements  revealed  that  the  spectral  features  characteristic  of 
screech  and  broadband  shock-associated  noise  are  present 
even  at  the  design  condition.2-  Large  Eddy  Simulation  (LES) 
by  the  present  authors  have  shown  the  presence  of  shocks 
inside  the  conical  C-D  nozzle  and  the  presence  of  shock  dia¬ 
monds  even  at  the  design  condition.-’ 

Narrow-band  spectra  measured  for  the  same  nozzle 
revealed  details  of  the  spectral  features  of  screech  and  broad¬ 
band  shock-associated  noise.24'25 

The  present  paper  expands  upon  the  initial  evidence  pre¬ 
sented  by  Munday  et  «/.22and  Liu  et  al.23  in  improving  the 
quality  of  the  flow  measurements  with  improved  shadow¬ 
graph  images  and  by  adding  Particle  Image  Velocimetry 
(PIV)  data.  The  present  paper  also  nearly  doubles  the  domain 
over  which  near-field  pressure  fluctuations  are  measured  and 
expands  the  number  of  operating  conditions  (Mj  or  NPR) 
over  which  far-field  acoustic  data  is  presented.  The  grid  on 
which  the  LES  has  been  performed  has  been  further  refined 
in  the  region  around  the  nozzle  exit  to  better  resolve  the 
shear  layer  and  shock  structures. 

B.  Experimental  facility  and  procedure 

The  three  conical  C-D  nozzles  employed  in  this  study 
are  simplified  static  models  of  actual  variable  geometry  noz¬ 
zles.  Rather  building  a  laboratory  model  with  variable 
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FIG.  1.  Nozzle  cross-section  with  schematic  of  a  double-diamond  shock 
cell  structure.  Dimensions  shown  are  in  milimeters.  The  three  nozzles 
examined  have  values  of  De  of  69.29,  72.85,  and  76.35  mm. 

geometry,  three  static  CCD  nozzles  have  been  produced  that 
approximate  the  real  engine  nozzles  in  three  different  config¬ 
urations  corresponding  to  three  different  operating  condi¬ 
tions.  The  laboratory  model  cross-sections  are  circular  rather 
than  polyhedral  so  the  surfaces  of  the  convergent  and  diver¬ 
gent  sections  are  conical.  There  are  no  hinges  at  the  throat, 
but  the  sharp  character  of  the  throat  is  retained.  The  geome¬ 
try  of  the  nozzles  is  shown  in  Figure  1 .  All  three  nozzles  are 
identical  upstream  of  the  throat.  They  all  have  a  throat  diam¬ 
eter,  Dt,  of  72.84  mm  and  expansion  sections  73.56  mm  in 
axial  length.  They  have  exit  diameters,  De,  of  69.29,  72.85, 
and  76.35  mm  giving  the  nozzles  area  ratios  of  1.067,  1.181, 
and  1.295  and  design  Mach  numbers,  Md,  of  1.3,  1.5,  and 
1.65,  respectively.  The  models  are  mounted  in  the  University 
of  Cincinnati  Acoustic  Test  Facility  (UC-ATF)  which  is  a 
7.32  x  7.62  x  3.35  m  test  chamber  which  has  been  acousti¬ 
cally  treated  to  be  anechoic  down  to  500  Hz.  This  chamber 
includes  fixtures  for  far-held  acoustics  and  PIV  flow 
measurements. 

C.  Experimental  and  numerical  test  conditions 

Three  experimental  nozzles  are  used  having  design 
Mach  numbers,  Md=1.3,  1.5,  and  1.65.  Each  nozzle  was 
operated  at  multiple  nozzle  pressure  ratios  (NPRs)  represent¬ 
ing  underexpanded,  design  condition,  and  overexpanded 
conditions.  All  cases  were  run  with  an  unheated  total  temper¬ 
ature  (T0  =  300  K). 

LES  results  presented  here  correspond  to  the  underex¬ 
panded  jet  condition  with  the  geometry  corresponding  to  the 
Md  =1.5  nozzle  with  total  temperature  T0  =  300  K. 

The  nozzle  pressure  ratios  used,  the  corresponding  fully 
expanded  jet  Mach  numbers,  Mj,  and  other  parameters  of  inter¬ 
est  are  shown  in  Table  I.  The  Reynolds  number,  Rcj,  is  based 
on  Uj  and  Dj,  where  Uj  is  the  jet  velocity  and  D,  is  the  fully 
expanded  jet  diameter.  /fd  is  a  shock  strength  parameter  defined 
as  ±  (IMj2  —  Md2l)0'5  where  the  plus  sign  is  taken  for  underex¬ 
panded  cases  when  Mj  is  greater  than  Md,  and  the  minus  sign  is 
taken  for  overexpanded  cases  when  Mj  is  less  than  Md.  This  pa¬ 
rameter  reduces  to  the  fl  introduced  by  Harper-Bourne  and  Fi¬ 
scher9  for  convergent  nozzles  where  Md  =  1 .0. 

D.  Instrumentation 

The  high  gradient  features  in  the  flow  are  visualized  by 
the  shadowgraph  technique.  An  Oriel  66056  arc  lamp  is  used 
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TABLE  I.  Test  conditions. 


Md 

NPR 

Mj 

T0[K] 

Uj  [m/s] 

Dj  [mm] 

Re,  [xlO6] 

ft 

1.3 

2.50 

1.22 

300 

372.6 

68.31 

2.6 

-0.44 

1.3 

2.77 

1.30 

300 

390.1 

69.24 

2.9 

0.00 

1.3 

3.25 

1.41 

300 

415.1 

71.07 

3.3 

0.56 

1.3 

3.69 

1.50 

300 

433.2 

72.80 

3.6 

0.76 

1.5 

2.50 

1.22 

300 

372.6 

68.31 

2.6 

-0.87 

1.5 

3.00 

1.36 

300 

402.9 

70.10 

3.1 

-0.64 

1.5 

3.50 

1.47 

300 

425.8 

72.05 

3.5 

-0.31 

1.5 

3.67 

1.50 

300 

432.5 

72.72 

3.6 

0.00 

1.5 

4.00 

1.56 

300 

444.0 

74.01 

3.9 

0.42 

1.5 

4.50 

1.64 

300 

458.8 

75.93 

4.2 

0.66 

1.5 

5.00 

1.71 

300 

471.3 

77.79 

4.6 

0.82 

1.65 

3.30 

1.43 

300 

417.4 

71.26 

3.3 

-0.83 

1.65 

4.59 

1.65 

300 

461.2 

76.27 

4.3 

0.00 

1.65 

6.30 

1.86 

300 

496.5 

82.32 

5.5 

0.86 

for  illumination.  A  pair  of  305  mm  parabolic  first-surface 
mirrors  with  a  1.83  m  focal  length  is  employed  to  collimate 
the  light  before  the  model  and  then  to  refocus  the  beam  after. 
The  image  is  captured  with  an  exposure  time  of  10  ps  using 
a  LaVision  Imager  Intense  cross-correlation  CCD  camera 
with  1376  x  1040  pixel  resolution  and  12-bit  intensity  reso¬ 
lution.  This  gives  a  spatial  resolution  on  the  order  of  0.004 
throat  diameters.  A  28-300  zoom  lens  is  mounted  to  the 
camera  allowing  optimization  of  the  field  of  view.  The  aper¬ 
ture  is  left  completely  open  and  exposure  is  controlled  by 
mounting  a  filter.  Averaging  100  images  eliminates  the  tur¬ 
bulence  and  gives  a  clear  view  of  the  shock  and  Prandtl- 
Meyer  waves. 

Detailed  flow-field  mapping  is  performed  by  PIV.  The 
PIV  system  is  built  by  LaVision  and  the  entire  PIV  suite,  laser 
and  cameras  are  mounted  on  a  traverse  which  allows  the  sys¬ 
tem  to  be  translated,  undisturbed,  to  any  streamwise  location 
allowing  many  fields  of  view  to  be  measured  without  the  loss 
of  time  in  changing  setups  and  without  the  uncertainties  which 
come  from  repeated  adjustment  of  components.  The  flow  is 
seeded  with  olive  oil  droplets  with  diameters  on  the  order  of 
1  /mi.  A  500  mJ  New  Wave  Research  Nd:YAG  double-pulse 
laser  is  passed  through  sheet-forming  optics  to  illuminate  the 
seed  and  the  images  are  captured  by  a  pair  of  LaVision  CCD 
cameras  with  1376  x  1040  pixel  resolution  and  12-bit  inten¬ 
sity  resolution.  Images  in  each  pair  are  separated  in  time  by 
1  ps,  giving  peak  particle  displacements  on  the  order  of  half  a 
millimeter.  The  image  scaling  is  7.09  pixels/mm  so  the  peak 
displacements  are  on  the  order  of  three  pixels.  The  image 
pairs  are  cross  correlated  using  a  multi-pass  method  beginning 
with  64  x  64  pixel  interrogation  windows  cascading  down  to 
4x4  pixel  windows  with  0%  overlap.  The  multi-pass  method 
applies  a  window  offset  for  each  subsequent  pass  based  on  the 
results  of  the  previous  pass,  improving  accuracy.  Assuming 
the  displacement  uncertainty  of  0.04  pixels  given  in  Ref.  26 
gives  an  uncertainty  in  velocity  of  5.6  m/s.  This  is  on  the  order 
of  0.0  lUj.  Individual  velocity  vectors  above  550  m/s  have 
been  discarded  as  spurious.  Discarded  vectors  are  not  replaced 
by  interpolation  in  the  final  pass.  Only  valid  vectors  are  con¬ 
sidered  for  averaging  or  other  statistical  analysis.  Histograms 
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FIG.  2.  Schematic  of  microphone  locations  for  near-field  and  far-field  meas¬ 
urements.  The  far-field  array  consists  of  10  microphones  ranging  from 
i jf  =  35°  to  150°,  each  located  3.8  m  from  the  center  of  the  nozzle  exit  plane. 
The  near-field  array  consists  of  a  rake  with  32  microphones  which  are  trav¬ 
ersed  through  13  locations  to  sweep  out  a  0.3  x  0.8  m  grid.  Two  such  grids  are 
swept  to  give  a  total  region  of  0.3  x  1 .6  m. 


of  velocities  in  the  individual  vector  fields  have  been  checked 
and  no  evidence  of  pixel  locking  was  found. 

Far-field  measurements  are  obtained  by  mounting  B&K 
Type  4954-B  quarter-inch  free-held  microphones  on  a  circular 
arc  56  throat  diameters  from  the  center  of  the  nozzle  exit 
plane  as  shown  in  Figure  2.  The  facility  is  described  in  detail 
by  Callender  et  air1  These  microphones  cover  an  arc  between 
the  angles  35°  and  150°  measured  from  the  upstream  direc¬ 
tion.  The  microphone  signals  are  low-pass  filtered  up  to  100 
kHz  and  recorded  at  204.8  kHz  for  10  s.  The  flow  variables 
are  monitored  at  4  Hz  in  order  to  establish  and  hold  the 
desired  operating  condition  of  the  model.  Variables  recorded 
include  the  ambient  conditions  in  the  chamber  as  well  as  the 
stagnation  temperature  and  pressure  of  the  jet. 

In  order  to  map  the  near-field  acoustics,  thirty-two 
microphones  are  mounted  at  2.54  cm  intervals  on  a  linear 
rake  which  is  mounted  on  a  large  three-axis  traverse  and  is 
then  traversed  radially  outward  by  steps  of  2.54  cm.  This 
produces  a  grid  of  31  x  12  measurement  locations.  Two 
31  x  12  grids  are  taken  to  extend  the  overall  region  to 
61  x  12  locations  covering  23  x  4.5  throat  diameters.  One 
column  of  measurement  locations  is  common  to  both  grids. 
The  location  of  the  measurement  grids  are  shown  in  Figure 
2.  The  rake  is  angled  at  10°  to  keep  the  microphones  out  of 
the  flow  as  the  jet  spreads.  The  resulting  grid  of  pressure 
fluctuations  can  then  be  analyzed  on  an  Overall  Sound 
Pressure  Level  (OASPL)  basis  or  for  particular  frequencies 
pertinent  to  the  different  mechanisms  and  can  be  mapped  to 
look  at  source  location  and  direction  of  propagation. 

Both  far-field  and  near-field  experimental  acoustic  data 
are  analyzed  in  blocks  of  4096  samples.  Each  block  is  win¬ 
dowed  with  a  Kaiser  filter  and  then  a  fast-Fourier-transform 
(FFT)  is  applied  yielding  narrow-band  spectra  with  a  fre¬ 
quency  resolution  of  50  Hz.  After  the  spectra  are  nondimen- 
sionalized  the  spacing  between  Strouhal  numbers  is 
approximately  0.125.  Finally  the  spectra  derived  from  all  the 
blocks  are  averaged.  OASPL  are  determined  from  the  time  do¬ 
main  data  by  applying  a  4th  order  high-pass  Butterworth  filter 
to  remove  content  below  500  Hz,  then  determining  the  Pms. 


E.  Numerical  approach 

The  unsteady  three-dimensional  inviscid  compressible 
flow  equations  are  solved  with  a  finite  element  code  FEFLO 
on  unstructured  tetrahedral  grids. 2S  This  code  is  capable  of 
accurately  representing  complicated  geometries,  such  as  the 
nozzle  geometries  used  in  this  work.  The  Flux-Corrected 
Transport  (FCT)  algorithm  is  used  for  the  spatial  discretiza¬ 
tion  and  a  fourth-order  Taylor-Galerkin  scheme  is  used  for 
the  time  integration.  FCT  is  ideal  for  the  shock  containing 
flows  because  it  is  high-order,  conservative,  monotone,  and 
positivity-preserving.24  No  explicit  subgrid  scale  model  is 
used  and  the  modeling  of  subgrid  scales  is  implicitly  pro¬ 
vided  by  the  embedded  flux  limiter.  The  present  simulations 
are  in  the  category  of  monotonically  integrated  large  eddy 
simulation  (MILES)  approach.30 

Characteristic  boundary  conditions  are  applied  to  both  far 
field  and  outflow  boundaries,  and  the  total  pressure  is  kept 
constant  at  the  inlet  of  the  nozzle.  A  slip  boundary  condition 
is  applied  to  the  nozzle  wall,  and  no  turbulence  is  added  at  the 
inflow  boundary.  The  cell  size  is  roughly  0.02  De  near  the 
nozzle  exit  and  it  gradually  increases  to  0.034De  in  the  wake 
region.  Since  the  jets  under  study  have  Reynolds  numbers 
greater  than  106,  the  shear  layer  near  the  nozzle  lip  is  thin  and 
a  fine -mesh  layer  with  a  cell  size  of  0.0035De  is  clustered  near 
the  nozzle  lip  and  up  to  0.5De  downstream  to  resolve  the  shear 
layer.  The  mesh  size  increases  progressively  and  becomes 
very  large  in  the  far-field  region.  A  buffer  region  is  added 
near  the  far-field  and  outflow  boundaries  to  avoid  wave  reflec¬ 
tions.  More  details  of  the  numerical  approach  and  its  valida¬ 
tion  can  be  found  in  Ref.  3 1 . 

After  the  flow  field  and  the  near-field  acoustic  fields 
reach  a  quasi-steady  state,  pressure  information  in  the  near 
field  are  collected  at  a  small  time  interval  (<  0.04  De/Uj), 
and  a  data  set  over  a  time  length  greater  than  900  De/Uj  is 
collected.  This  data  set  is  divided  into  three  blocks  and  each 
one  has  8192  samples  with  a  frequency  resolution  of  20  Hz. 
The  final  acoustic  data  are  calculated  by  ensemble  averaging 
data  over  all  blocks.  Both  OASPL  and  the  space-time  corre¬ 
lations  are  evaluated  over  the  entire  data  set  and  the  near¬ 
field  OASPL  are  evaluated  from  Prms 

II.  RESULTS  AND  DISCUSSION 
A.  Flow  structure 

LES  reveals  the  flow  inside  the  nozzle  where  it  is  difficult 
to  make  any  unobtrusive  flow  measurements.  Figure  3  shows 
mean  Mach  number  contours  inside  the  nozzle  and  in  the 
region  just  past  the  nozzle  exit  for  Md=  1.5  and  Mj  =  1.56. 
The  Mach  contours  reveal  that  the  sharp  throat  produces  an 
oblique  shock  inside  the  nozzle  and  forms  a  Mach  disk  near 
the  centerline.  For  this  case,  the  Mach  disk  lies  at 
x/De=— 0.35  which  is  inside  the  nozzle  and  the  interface 
between  the  oblique  shocks  and  the  Mach  disk  casts  a  slip  line 
near  the  jet  centerline  with  a  diameter  of  r/De  =  0.06.  The 
onward-running  oblique  shock  extends  from  this  internal 
Mach  disk,  past  the  nozzle  exit  and  passes  through  the  oblique 
shock  anchored  at  the  nozzle  lip.  It  reflects  from  the  shear 
layer  outside  the  nozzle  downstream  from  the  nozzle  lip  at 
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FIG.  3.  Mean  Mach  numbers  from  LES  inside  the  nozzle  and  near  the  exit 
for  Md  =  1.50,  Mj  =  1.56. 


x/De  =  0.07.  The  wave  reflected  from  this  throat  shock  is  an 
expansion  wave  that  runs  parallel  to  the  lip  wave  towards  the 
centerline. 

The  case  shown  in  Figure  3  (Md  =  1.5,  Mj=  1.56)  is  a 
mildly  underexpanded  case.  The  classical  theory  would  sug¬ 
gest  that  overpressure  will  cause  the  flow  to  turn  outward  in 
an  expansion  fan  which  would  be  anchored  at  the  nozzle  lip 
and  run  inward  toward  the  axis.32  Since  this  nozzle  is  coni¬ 
cal,  the  flow  exits  with  a  radially  outward  component  of  ve¬ 
locity  and  must  turn  inward,  despite  the  underexpanded 
condition,  so  the  wave  running  inward  from  the  lip  is  a 
shock.  This  can  be  seen  by  the  drop  in  Mach  number  from 
the  upstream  side  of  the  lip  wave  to  the  downstream  side. 

The  flow  pattern  shown  in  Figure  3  differs  from  a  classi¬ 
cal  supersonic  jet  in  two  ways,  both  due  to  the  particular 
details  of  the  nozzle.  First,  the  conical  expansion  forces  an 
inward  turning  of  the  flow  near  the  nozzle  lip  for  all  cases 
near  the  design  condition.  The  flow  leaves  the  nozzle  with  an 
outward  radial  component  of  velocity  and  this  must  turn 
inward.  Secondly,  the  sharp  throat  generates  an  oblique  shock 
regardless  of  the  condition  outside  the  nozzle.  The  combined 
result  of  these  two  factors  is  that  there  are  two  overlapping 


sets  of  shock  diamonds  produced  and,  as  will  be  shown,  these 
double  diamonds  are  present  for  the  design  condition  as  well 
as  for  underexpanded  and  overexpanded  conditions.  This 
double-diamond  feature  is  depicted  schematically  in  Figure  1 . 

In  Figure  4  shadowgraph  and  PIV  results  are  compared. 
The  case  is  Md  =  1.5  and  Mj  =  1.56.  Instantaneous  and  aver¬ 
aged  shadowgraphs  are  shown  in  Figs.  4(a)  and  4(b),  respec¬ 
tively.  Instantaneous  and  averaged  axial  velocities  from  PIV 
are  shown  in  Figs.  4(c)  and  4(d). 

Figure  4(a)  shows  an  instantaneous  shadowgraph  image 
of  the  jet  taken  with  an  exposure  time  of  10  ps.  The  turbu¬ 
lence  near  the  nozzle  appears  fine  grained,  but  the  scale  of 
the  turbulence  increases  in  the  downstream  direction.  The 
outer  edges  of  the  jet  are  ragged  with  large  scale  structures 
and  as  these  convect  they  produce  Mach  wave  radiation 
which  can  be  seen  propagating  at  a  downstream  angle  to 
the  main  jet.  Oblique  shocks  can  be  seen  emerging  from 
the  nozzle  exit  and  reflecting  from  the  shear  layer  at 
x/De  =  T0/De  =  0.07,  where  T0  is  the  axial  position  of  the 
reflected  throat  wave  as  shown  in  Figure  1 . 

These  reflections  make  bright  spots  in  the  shadowgraph 
at  the  upper  and  lower  shear  layers  and  a  bright  vertical  line 
connecting  these  spots.  This  is  due  to  the  integration  along 
the  line  of  sight.  The  jet  geometry  is  axisymmetric  and  so 
are  the  oblique  shocks  and  their  reflections.  The  bright  spots 
at  the  top  and  bottom  of  the  image  are  reflections  viewed  in 
profile.  Where  the  bright  line  crosses  the  jet  centerline  the 
camera  is  viewing  the  reflection  through  the  near  and  far 
edges  of  the  jet.  The  reflection  of  the  shock  is  outside  the 
field  of  view  of  the  PIV  images,  but  can  be  compared  with 
LES  Mach  number  contours  in  Figure  3.  The  LES  image  is  a 
section  in  an  axial  plane  and  is  not  line  integrated.  The 
reflection  of  the  throat  shock  is  clearly  visible  in  Figure  3  but 
the  vertical  line  is  not  present  since  its  origin  is  not  in  the 
plane  computed  by  the  LES. 

The  first  complete  double  diamond  can  be  seen  in  Figure 
4(a),  but  as  the  turbulence  is  more  pronounced  farther  down¬ 
stream,  the  second  cell  is  largely  obscured.  In  order  to  reveal 
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FIG.  4.  Flow  field  representations  by 
shadowgraph  and  PIV;  Mj  =  1.56 
Mj  —  1 .50.  (a)  Instantaneous  shadow¬ 
graph.  (b)  Average  of  100  shadowgraph 
images,  (c)  Instantaneous  PIV.  (d)  Aver¬ 
age  of  500  PIV  image  pairs.  PIV  veloc¬ 
ities  below  0.1  u/Uj  are  not  shown. 
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the  relatively  stationary  shock  cell  structures,  100  shadow¬ 
graph  images  were  averaged.  This  process  averages  out  the 
turbulence  and  leaves  the  shock  structures  visible.  The  aver¬ 
aged  shadowgraph  image  is  shown  in  Figure  4(b).  In  addition 
to  the  throat  shock  emerging  from  the  nozzle  and  the  first 
cell  of  the  double  diamond,  the  averaging  brings  out  the 
reflection  of  the  lip  shock  at  Li/De=  1.1  and  the  reflection  of 
the  throat  shock  at  T ]/De  =  1 .4,  and  the  next  set  of  reflections 
at  L2/De  =  2.4  and  T2/De  =  2.8.  Ln  and  Tn  are  the  axial  coor¬ 
dinates  of  the  nth  lip  and  throat  wave  reflections,  respec¬ 
tively  (see  Figure  1).  The  shock  structures  become  less 
sharply  imaged  as  they  get  farther  from  the  nozzle.  This  is 
due  to  the  fact  that  they  are  not  completely  stationary.  The 
shear  layer  is  in  motion  as  large-scale  structures  pass  along 
it.  Each  reflection  of  a  wave  off  the  shear  layer  is  reflecting 
off  a  moving  entity  so  the  reflected  waves  are  moved  by  the 
reflections.  This  effect  is  cumulative  so  as  the  number  of 
reflections  increases,  the  degree  of  motion  of  the  shocks  and 
expansion  waves  increases,  and  the  averaged  images  of  these 
moving  waves  become  increasingly  smeared  with  each  wave 
reflection. 

Another  feature  revealed  by  the  averaging  of  the  shad¬ 
owgraph  images  in  Figure  4(b)  is  the  slip  line  cast  from  the 
Mach  disk  inside  the  nozzle.  This  appears  as  a  dark  horizon¬ 
tal  feature  running  along  the  centerline  of  the  jet  images.  The 
width  of  the  slip  line,  w,  is  the  same  as  the  diameter  of  the 
internal  Mach  disk.  The  presence  of  the  slip  line  in  the  shad¬ 
owgraph,  and  PIV  to  follow,  provides  confirmation  of  the 
presence  of  the  Mach  disk  predicted  by  the  LES  in  Figure  3. 

The  bright  vertical  line  associated  with  the  throat  wave 
reflection  at  T0  observed  in  the  instantaneous  shadowgraph 
image  is  present  in  the  averaged  image  as  well.  The  averag¬ 
ing  out  of  the  turbulence  reveals  similar,  though  less  sharp 
vertical  lines  associated  with  the  reflections  at  L1;  Th  L2, 
and  T2  also.  As  previously  discussed,  these  are  artifacts  of 
the  fact  that  the  shadowgraph  image  is  a  line  integrating 
result.  Comparing  these  features  in  the  averaged  shadow¬ 
graph  in  Fig.  4(b)  with  the  corresponding  averaged  PIV  in 
Fig.  4(d)  shows  that  while  the  reflections  are  present,  the  ver¬ 
tical  lines  are  not. 

PIV  in  an  axial  plane  is  shown  in  Figures  4(c)  and  4(d). 
A  single  instantaneous  measurement  from  a  single  image 
pair  is  shown  in  Fig.  4(c)  and  an  average  of  500  such  instan¬ 
taneous  measurements  is  shown  in  Fig.  4(d).  While  PIV  has 
the  benefit,  relative  to  shadowgraph,  of  not  being  line  inte¬ 
grating,  it  has  a  disadvantage  in  regions  near  a  solid  surface 
like  the  nozzle  because  light  reflected  from  the  hardware 
over-saturates  the  camera.  For  this  reason  the  PIV  domain 
begins  near  x/De  =  0.25,  away  from  the  reflecting  hardware. 
The  plots  shown  are  of  the  axial  component  of  velocity  and 
velocities  below  u/Uj  =  0.1  have  been  omitted  in  order  to 
mask  the  background  noise. 

The  instantaneous  measurement  in  c  shows  that  the 
larger  scale  turbulence  observed  in  the  instantaneous  shad¬ 
owgraph  is  confined  to  the  shear  layer.  In  the  shadowgraph 
one  can  see  large  scale  turbulence  over  the  entire  jet  because 
of  line  integration.  PIV  shows  a  cut  through  an  axial  plane 
and  shows  that  the  variation  within  the  jet  is  small.  There  is 
large  variation,  however,  in  the  shear  layers.  These  are  thin 


at  the  nozzle  exit  and  thicken  significantly  in  the  downstream 
direction.  Large  variation  in  velocity  is  apparent  in  the  thick¬ 
ened  shear  layer  as  high-speed  eddies  entrain  relatively  still 
air  from  the  jet’s  surroundings.  Pockets  of  velocity  below 
u/Uj  =  0.1  (shown  in  white)  are  drawn  deep  into  the  shear 
layer. 

The  velocity  differences  across  the  shocks  and  expan¬ 
sion  waves  in  the  shock  cells  are  small  enough  to  be  compa¬ 
rable  to  the  random  variation  in  velocity  so  the  shock  cells 
are  not  very  pronounced  in  the  instantaneous  velocity  plots. 
Averaging  500  measurements  suppresses  this  variation  and 
the  mean  velocity  plot  in  Figure  4(d)  clearly  shows  the 
double-diamond  structure  and  the  presence  of  the  slip  lines 
around  the  centerline.  The  initial  reflection  of  the  throat 
shock  occurs  outside  the  measurement  domain,  but  all  the 
remaining  features  previously  observed  are  present. 

As  with  the  averaged  shadowgraph,  the  first  shocks  are 
the  sharpest  in  character.  Later  shocks  become  somewhat 
smeared  since  they  fluctuate  in  location  due  to  accumulated 
variation  from  additional  reflections  off  the  shear  layer.  The 
locations  of  the  later  waves  should  be  understood  to  be  aver¬ 
age  locations  since  the  waves  are  in  motion. 

The  wake  of  the  internal  Mach  disk,  bounded  by  two 
slip  lines  is  evident  throughout  the  measurement  domain. 
The  width  of  this  slip  line,  w,  was  measured  for  all  three  noz¬ 
zles  and  over  a  range  of  jet  Mach  numbers  and  it  was  found 
that  it  is  chiefly  a  function  of  design  Mach  number.  Increased 
Md  and  increased  area  ratio  lead  to  decreasing  slip  line 
width,  and  thus  decreasing  diameter  of  the  internal  Mach 
disk.  The  slip  line  width  w/De  is  0.09,  0.06,  and  0.045  for  Md 
of  1.3,  1.5,  and  1.65,  respectively.  There  is  also  a  very  weak 
dependence  on  fully  expanded  jet  Mach  number,  showing 
reduced  w/De  with  increased  Mj. 

Comparison  of  Figure  4(d)  with  Figure  3  shows  that  the 
features  associated  with  the  double-diamond  structure  pre¬ 
dicted  by  LES  agree  well  with  the  PIV  and  shadowgraph. 
The  shock  and  Prandtl-Meyer  wave  locations  and  the  posi¬ 
tion  of  the  wake  from  the  internal  Mach  disk  match  and  the 
LES  provides  a  much  clearer  image  of  these  than  the  PIV. 

The  axial  velocity  comparison  between  PIV  and  LES  is 
shown  in  Fig.  5,  where  the  distributions  along  the  centerline 
and  along  a  line  at  r/De  =  0.25  are  presented.  The  shock-cell 
shape  and  location  are  in  a  very  good  agreement  between 
these  two  predictions,  and  the  difference  of  the  shock-cell 
locations  is  within  0.04De.  The  major  difference  between 
these  two  predictions  is  the  magnitude  of  the  velocity;  PIV 
measurement  is  roughly  4.0%  higher  than  that  of  the  LES 
prediction  away  from  the  jet  center. 

Figure  6  shows  shadowgraph  images  for  a  range  of  NPRs 
including  two  overexpanded  conditions,  the  design  condition 
and  two  underexpanded  conditions.  There  are  double¬ 
diamonds  apparent  in  all  the  shadowgraph  images  including 
the  design  condition  (Mj  =  Md=1.5).  This  shows  that  the 
CCD  nozzles  studied  in  this  paper  do  not  achieve  a  shock-free 
condition  at  the  design  condition.  A  shadowgraph  movie  was 
taken  of  the  nozzle  as  the  set-point  was  slowly  varied  through 
a  range  from  Mj  =  1.22  through  Mj=  1.71  (NPR=2.0  up  to 
NPR=5.0),  including  the  points  around  Mj  =1.5  where 
quasi- 1-D  theory  predicts  perfect  expansion.  Careful 
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FIG.  5.  PIV  and  LES  velocities.  Mi  =1.56  Md=1.50.  (a)  r/De  =  0. 
(b)  r/De  =  0.25. 


observation  of  the  movie  reveals  that  there  is  no  condition  in 
this  range  of  Mj  in  which  the  jet  is  shock  free. 

For  underexpanded  conditions  a  decrease  in  Mj  shortens 
the  shock  cells  and  also  moves  the  two  diamond  structures 
closer  to  one  another  (L,  closer  to  T,)  until  at  the  lowest  con¬ 
dition  the  two  sets  of  waves  appear  to  have  coalesced  by  the 
first  lip  wave  reflection.  By  Mj=1.22  the  lip  shock  has 
formed  a  Mach  disk  three  times  the  diameter  of  the  internal 
Mach  disk. 

For  overexpanded  cases,  increased  Mj  causes  the  reflec¬ 
tions  to  grow  closer  to  one  another  as  the  shock  cells  elon¬ 
gate.  As  the  initial  Prandtl-Meyer  fan  from  the  lip 
strengthens  it  becomes  physically  thicker  and  comes  to  en¬ 
velop  the  throat  waves  entirely  for  the  highest  values  of  Mj. 

One  feature  which  is  highly  sensitive  to  Md  but  rela¬ 
tively  insensitive  to  Mj  is  the  position  of  the  initial  throat 
wave  reflection  T0.  For  Md  =  1.5  it  is  0.06  or  0.07  De  for  all 
cases  except  for  Mj  =  1.22  where  To  is  slightly  closer  to  the 
exit  (0.04De).  For  the  larger  nozzle  Md=  1.65  the  position  of 
the  initial  throat  reflection  is  farther  from  the  nozzle  at 
To/De  =  0.22,  0.29,  and  0.32  for  Mj  of  1.43,  1.65,  and  1.86, 
respectively.  The  smaller  nozzle,  with  Md=  1.3,  has  the  ini¬ 
tial  reflection  of  the  throat  wave  inside  the  nozzle.  Oil  depo¬ 
sition  within  the  nozzle  observed  after  running  reveals  the 
location  of  the  reflection  to  be  at  T0/De  =—  0.14. 

Figure  7  shows  instantaneous  static  pressure  distribu¬ 
tions  predicted  by  LES  for  a  mildly  underexpanded  case 
(Md=  1.5,  Mj  =  1 .56).  The  high  and  low  pressure  regions 
alternating  along  the  body  of  the  jet  mark  the  shock  cells. 
Initially  these  clearly  have  an  overall  diamond  shape. 
Beyond  x/De  =  4  or  5  the  instantaneous  images  are  less  dis¬ 
tinctly  shaped  and  they  begin  to  be  shifted  away  from  the  jet 
centerline  as  the  jet  oscillates.  By  x/De=  10  or  12  the  clear 
alternation  of  high  and  low  pressures  in  the  axial  direction 


FIG.  6.  (Color  online)  Convergent-divergent  nozzle  (Md=1.5)  shadow¬ 
graph.  One  hundred  images  averaged  to  average  out  turbulence  and  empha¬ 
size  shock  structures. 


ends.  This  marks  the  end  of  the  potential  core.  Beyond  this 
point  the  mixing  has  penetrated  all  the  way  to  the  jet  center 
and  the  mixing  process  loses  its  cylindrical  character. 

Near  the  nozzle  (x/De  =  0  to  3)  small  pressure  pockets 
can  be  seen  along  the  edge  of  the  jet.  These  are  signatures  of 
large-scale  structures  originating  at  the  nozzle  lip  and  con- 
vecting  along  the  shear  layer.  They  interact  with  the  shock 
cells  and  grow  larger  as  they  move  downstream.  Beyond 
x/De  =  3  they  become  comparable  in  size  to  the  pressure 
pockets  in  the  body  of  the  jet  caused  by  the  shock  cells. 
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Pressure  waves  with  large  intensities  radiate  from  the 
region  where  the  convecting  large-scale  structures  have 
grown  (x/De  >  3)  where  the  source  is  the  vigorous  interaction 
of  the  large-scale  structures  with  the  shock  cells,  producing 
screech  and  broadband  shock-associated  noises.  Large  pres¬ 
sure  waves  also  radiate  from  the  mixing  region  (x/De  >  10) 
where  mixing  noise  is  the  dominant  feature. 

Axial  locations  (L,  and  T,)  of  the  reflections  of  both  the 
lip  waves  and  tip  waves  have  been  extracted  from  the  PIV 
measurements.  The  throat  and  lip  reflections  are  only  distin¬ 
guishable  for  the  first  two  shock  cells.  Beyond  this  point  the 
smearing  due  to  shock  motion  blends  the  two.  The  shock  cell 
spacing  is  nearly  uniform  over  the  first  five  cells.  Beyond 
this  the  cells  become  shorter  due  to  the  thickening  of 
the  shear  layer.  This  thickening  reduces  the  diameter  of  the 
potential  region  and  reduces  the  sonic  diameter  at  which  the 
reflections  occur.  Reducing  the  diameter  of  the  cells  also 
reduces  their  length.  In  order  to  give  a  single  measure  Ls  of 
shock  cell  length,  an  average  has  been  taken  over  the  first 
five  shock  cells. 

The  shock  cell  spacing,  Ls  from  PIV  and  LES  is  shown 
in  Figure  8.  Also  shown  are  the  theoretical  values  using  the 
Prandtl-Pack  relation,33 

Ls  =  n(Mj  -  l)l/2Dj/ni,  where  /q  =  2.40483 

as  well  as  published  data  for  a  range  of  design  Mach  num¬ 
bers  from  1 .0  to  2.0.  Nearly  all  of  the  C-D  nozzles  previously 
published  have  been  of  the  smoothly  varying  type,  except  in 
Seiner  and  Norum20  and  in  Kuo,  Veltin  and  McLaughlin24  in 
which  military  style  C-D  nozzles  were  investigated.  The 
Prandtl-Pack  relation  is  generally  a  good  fit  for  convergent 
nozzles  (Md=1.0)  depicted  by  small  x’s.  But  it  generally 
over-predicts  Ls  for  higher  values  of  Md.  This  over- 


»  Ref  20,  M„  =  1.0 

■  Ref  1 1 ,  Ma  =  1 .0 

o  PIV.  M0  =  1.3 

Shadowgraph,  Mri  =  1.5 
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»  C/L  pressure,  Md  =  1.5 

<  LES,  Md  =  1.5 

-  Ref  34  conical,  M0  =  1.5 

t  Ref  34  smooth,  M0  =  1.5 
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*  PIV,  Md  =  1.65 

,  Ref  20.  Md  =  2,0 

r  Ref  11,  Md  =  2.0 

-  Prandtl-Pack,  Ref  32 

FIG.  8.  Shock  Cell  size,  Ls  as  a  function  of  Mj. 


prediction  becomes  worse  as  Mj  increases.  The  general  trend 
is  that  Ls  increases  with  Mj,  but  for  a  given  value  of  Mj  the 
shock  cell  size  decreases  with  Md.  The  values  of  Ls  extracted 
from  PIV,  from  Shadowgraph,  from  centerline  pressure 
measurements,  and  from  LES  all  agree  fairly  well  with  one 
another  and  agree  with  other  published  results  for  nozzles 
with  design  Mach  numbers  near  1.5.  It  can  be  seen  from  this 
that  there  is  not  any  significant  difference  in  Ls  introduced 
by  the  use  of  the  more  realistic  CCD  nozzles. 


B.  Acoustic  emissions 

Far-field  acoustic  spectra  for  the  Md=1.5  nozzle  are 
shown  in  Figure  9  for  three  conditions.  M,  =  1.47,  1.50,  and 
1.56  representing  mildly  overexpanded,  design  condition, 
and  underexpanded.  The  frequencies  are  expressed  in  terms 
of  a  Strouhal  number  based  on  Uj  and  Dj.  All  three  show  sig¬ 
natures  of  shock  associated  noise.  The  narrow  peak  between 
St  =  0.2  and  0.3  which  radiates  most  strongly  in  the  upstream 
direction  {ip  <  90°)  is  a  screech  peak  and  is  as  prominent  at 
the  design  condition  as  the  two  off-design  cases.  All  three 
sets  also  show  the  characteristic  signature  of  broad-band 
shock-associated  noise  (BBSN).  The  BBSN  peak  is  much 
broader  than  the  screech  peak,  is  centered  at  St  =  0.3  or  0.4 
at  t/f  =  35°,  and  shifts  to  higher  Strouhal  numbers,  grows 
broader,  and  lower  as  it  shifts  to  more  aft  observers.  By 
ip  =  140°  the  BBSN  peak  has  shifted  to  St  =  1.1  and  is  barely 
perceptible  above  the  background  mixing  noise.  At  i/f  =  150° 
the  BBSN  peak  is  negligible. 

This  similarity  between  the  spectra  of  the  design  condi¬ 
tion  to  those  of  off  design  conditions  is  characteristic  of  the 
conical  C-D  nozzles  used  in  this  study.  Traditional  smoothly 
varying  C-D  nozzles  become  shock  free  at  their  design  con¬ 
dition,  and  therefore  do  not  produce  screech  or  BBSN  at  or 
near  design  condition.  Operating  a  jet  on-design  can  there¬ 
fore  offer  significant  reduction  in  noise  for  smoothly  varying 
nozzles.  For  the  practical  engine  nozzles  used  on  high- 
performance  jet  aircraft  this  saving  is  not  available  even 
when  the  engine  area  ratio  is  closely  matched  to  the  operat¬ 
ing  condition. 

Spectra  for  a  far-field  microphone  at  t p  =  35°  are  shown 
for  all  three  nozzles  as  contour  plots  in  Figure  10.  The  design 
condition  for  each  nozzle  is  indicated  by  the  horizontal  black 
line.  These  plots  were  made  by  sweeping  slowly  through  val¬ 
ues  of  Mj  so  that  any  shock  free  spectra  could  be  evident.  No 
condition  was  found,  however,  for  which  there  was  no  signa¬ 
ture  of  shock  noise.  The  dominant  ridge  in  the  plot  is  the 
screech  tone.  All  three  nozzles  showed  strong  screech  away 
from  the  design  condition  with  a  diminution  in  intensity  and 
a  frequency  shift  at  the  point  at  which  the  screech  switched 
modes.  The  Md=1.3  nozzle  shows  two  such  shifts  at 
Mj=  1.25  and  1.4.  The  Md=  1.5  and  1.65  nozzles  only  ex¬ 
hibit  one  such  shift  near  the  design  condition.  The  broader 
and  lower  ridge  to  the  right  of  the  screech  tone  is  the  BBSN 
peak  for  each  nozzle.  This  BBSN  ridge  is  significant  near  the 
design  condition  for  all  three  nozzles.  There  is  a  slight  dimi¬ 
nution  in  intensity  for  the  Md=1.3  nozzle  near  Mj=1.4 
where  the  screech  intensity  is  also  at  a  minimum,  but  for  the 
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FIG.  9.  Far-field  acoustic  measurements  for  the  conical  convergent- 
divergent  nozzle  with  Md  =  1.5.  Each  trace  is  for  a  different  angle,  i//,  meas¬ 
ured  from  the  upstream  direction.  All  measurements  were  taken  at  a  distance 
of  56  throat  diameters,  (a)  Mj=  1.47,  overexpanded,  (b)  Mj  =  1.5,  perfectly 
expanded,  and  (c)  Mj  =  1.56,  underexpanded. 

M(l  =  1 .5  and  1.65  nozzles  there  is  no  reduction  in  BBSN 
whatsoever. 

The  peak  frequency  of  the  broadband  shock  associated 
noise  is  shown  in  Figure  1 1  for  an  observer  at  i ]/  =  90°.  These 
values  are  compared  with  published  results  for  smoothly 
varying  nozzles  and  for  the  military  C-D  nozzles  tested  by 
McLaughlen’s  group.34  The  data  for  all  nozzles  collapse, 
whether  smooth  or  conical,  for  all  design  Mach  numbers. 
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FIG.  10.  Far-field  spectra  at  1^  =  35°.  (a)  Md=1.3.  (b)  Md=1.5. 
(c)  Md  =  1.65. 

There  are  also  two  theoretical  curves  due  to  Tam  and 
Tanna.10  Equation  of  Tam  and  Tanna  is  based  on  a  different 
derivation  than  Harper-Bourne  and  Fisher’s4  but  the  equation 
for  the  principal  BBSN  peak  is  the  same. 


Ls(  1  H - -cos  i/t) 

@00 

where  Uc  the  convective  velocity  is  taken  to  be  0.7  Ur 
Tam’s  derivation  depends  on  knowledge  of  the  shock  cell 
spacing,  Ls,  for  which  he  uses  Prandtl-Pack’s  equation 
above.  The  Tam-Prandtl-Pack  line  on  the  plot  uses  Prandtl- 
Pack  for  this.  Tarn-Shadowgraph  uses  the  shock  cell  spacing 
read  from  the  shadowgraph  images  in  the  present  study.  The 
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Tam-Prandtl-Pack 
- Tam -Shadowgraph 

FIG.  11.  Peak  BBSN  frequencies  with  varying  Mj  for  an  observer  \jj  =  90°. 

Tam-Prandtl-Pack  equation  under-predicts  the  BBSN  fre¬ 
quency  for  nearly  all  cases,  but  correctly  represents  the 
trends.  Substituting  a  measured  value  for  Ls  improves  the 
Tam  and  Tanna  equation  considerably,  but  it  remains  on  the 
low  side.  The  trend  in  present  measurements  is  for  the  BBSN 
frequency  to  rise  with  Md  especially  at  lower  values  of  Mr 
Figure  12  shows  a  similar  plot  of  screech  frequency 
measurements,  compared  to  published  data  and  to  theoretical 
curves  expanded  to  predict  screech.  The  agreement  among 
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FIG.  12.  Strouhal  number  of  screech  peak  as  a  function  of  Mj. 
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FIG.  13.  Variation  in  OASPL  with  fiA  =  ±  (IMj2  -  Md2l)a5. 

different  investigators  is  good,  and  the  trend  is  for  increasing 
Screech  frequency  with  increasing  design  Mach  number 
especially  for  lower  values  of  Mj.  Tam,  Seiner,  and  Yu 11 
have  proposed  an  explanation  of  the  relationship  between 
screech  and  BBSN  which  holds  that  the  upstream  propagat¬ 
ing  component  of  the  screech  feed-back  loop  is  in  fact  the 
BBSN  noise,  and  this  the  screech  frequency  will  be  the 
BBSN  frequency  as  propagated  to  i/r  =  O' .  As  for  the  BBSN 
frequency,  so  the  screech  frequency  is  under-predicted  by 
Tam-Prandtl-Pack,  but  is  better  with  Tarn-Shadowgraph. 
Screech  is  predicted  by  our  LES  and  the  agreement  with 
experiment  is  good. 

The  CCD  nozzles  under  study  agree  with  smoothly  vary¬ 
ing  nozzles  in  the  frequencies  of  screech  and  BBSN,  but  the 
key  difference  acoustically  is  the  presence  of  shock-associated 
noise  near  the  design  condition  for  conical  C-D  nozzles  where 
it  is  absent  for  smoothly  varying  nozzles.  Plots  of  OASPL  as  a 
function  of  /fd  (Figure  13)  show  that  there  is  a  broad  region  of 
diminished  OASPL  around  the  design  condition,  but  as  dis¬ 
cussed  in  Fig.  9  this  is  due  to  reduction  in  screech  and  there  is 
no  reduction  in  BBSN.  Other  published  data  co-plotted  here 
show  that  for  traditional  C-D  nozzles  the  reduction  in  OASPL 
near  the  design  condition  is  more  significant  than  the  present 
study.  Seiner  and  Norum-0  data  of  the  “F15  nozzle”  which  is 
likely  similar  to  the  conical  C-D  nozzles  shows  similar  trends 
to  the  present  results. 


C.  Acoustic  source  locations 

Contour  plots  of  OASPL  and  Sound  Pressure  Level 
(SPL)  for  select  Strouhal  numbers  are  presented  in  Figure  14 
for  Md=1.5  and  Mj  =  1 .56.  Both  measurement  and  LES 
results  are  shown.  OASPL  and  SPL  from  LES  occupy  the  rec¬ 
tangular  region  from  r/De  =  0  to  3  and  x/De  =  0  to  25  except 
for  the  region  covered  by  near-held  microphone  measure¬ 
ments.  The  microphone  data  are  shown  in  the  rectangular 
region  angled  10°  from  the  jet  axis  and  extending  to  r/De  =  9 


116102-11 


Flow  structure  and  acoustics  of  supersonic  jets 


Phys.  Fluids  23,  1 1 61 02  (201 1 ) 


FIG.  14.  Near-field  contour  maps  of 
OASPL  and  50  Hz  bands  for  a  slightly 
underexpanded  case,  Md=1.50, 
M|=  1.56.  Levels  are  in  dB.  (a)  Over  all 
sound  pressure  level,  (b)  Mixing  noise, 
St  =  0.166.  (c)  Screech,  St  =  0.266.  (d) 
BBSN,  St  =  0.566. 
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and  x/De  =  22.  When  the  two  regions  overlap,  near-field  mea¬ 
surement  is  shown. 

Figure  14(a)  shows  the  OASPL.  First,  note  that  the  tran¬ 
sition  from  LES  to  measurement  is  very  smooth,  except  a 
very  small  region  near  the  nozzle  exit  where  the  inflow  tur¬ 
bulence  dominates  the  pressure  wave  radiation.  The  jet  is 
making  noise  across  the  entire  length  of  the  shock  train.  The 
maximum  OASPL  has  an  apparent  source  around  x/De=  10, 
but  high  levels  of  OASPL  range  from  x/De  =  2.5  to  13. 
Sound  pressure  levels  (SPLs)  in  50FIz  bands  around  selected 


frequencies  are  shown  in  the  remaining  subfigures.  Figure 
14(b)  shows  a  Strouhal  number  of  0.166.  This  is  in  a  fre¬ 
quency  range  below  which  no  shock  associated  noise  is  gen¬ 
erated,  so  it  represents  the  mixing  noise.  Again  a  very  good 
agreement  is  found  between  computations  and  measurement. 
Mixing  noise  appears  to  be  radiating  from  the  latter  part  of 
the  shock  train  and  extending  beyond  its  end.  Decibel  levels 
above  130  appear  ranging  from  x/De  =  7.5  to  17.5.  This 
apparent  source  region  includes  the  region  of  intense  mixing 
at  the  end  of  the  jet’s  potential  core.  Figure  14(c)  at  a 
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FIG.  15.  Space-time  fluctuating  pressure 
correlations  (p'(0, 2De,  To);  p'{x,  0.5/.),,. 

+  T))//7^2between  a  reference  point 
in  the  nozzle  exit  plane  at  two  diameters 
and  points  along  the  lip  line 
(r  =  time  ■  Dg/a^). 


Strouhal  number  of  0.266  is  the  screech  frequency  observed 
in  the  far- field.  This  plot  shows  the  signature  of  screech  in 
the  form  of  a  lobed  source  distribution  with  lobes  centered 
on  the  shock  reflections.  The  predictions  of  the  lobe  loca¬ 
tions  and  the  intensity  agree  with  measurement  well.  Figure 
14(d)  at  Strouhal  number  of  0.566  is  based  on  the  frequency 
of  the  BBSN  peak  observed  in  the  far-held  by  an  observer  at 
i//  =  90°.  This  distribution  shows  a  single  lobe  centered  on 
x/De=  10.5.  This  is  the  region  where  the  large  scale  struc¬ 
tures  encounter  the  latter  shock  cells.  For  this  frequency,  the 
overall  agreement  between  the  measurement  and  LES  pre¬ 
dictions  is  good,  but  it  appears  LES  under-predicts  the  inten¬ 
sity  slightly. 

Two-point  correlations  between  fluctuating  pressures  in 
the  near-  and  far-helds  have  been  used  to  reveal  source  loca¬ 
tions.3^  Since  the  LES  in  this  case  does  not  extend  to  the  far- 
held,  a  reference  point  is  chosen  in  the  near-held.  The  refer¬ 
ence  point  chosen  is  in  the  plane  of  the  nozzle  exit 
(x/De  =  0)  and  well  outside  the  jet  (r/De  =  2).  This  reference 
point  is  in  the  region  in  which  the  shock-associated  compo¬ 
nents  of  jet  noise  are  known  to  propagate.  Figure  15  shows 
correlations  of  instantaneous  pressure  between  the  reference 
point  and  points  along  the  lip  line  (r/De  =  0.5)  for  a  range  of 
axial  positions.  The  temporal  delay  t  ( time  ■  De/a^_)  is  on 
the  ordinate  and  x/De  is  on  the  abscissa.  The  peak  correlation 
regions  appear  as  streak  lines  with  alternate  positive  and  neg¬ 
ative  values  in  time,  indicating  that  pressure  fluctuations 
along  the  lip  line  with  large  correlation  values  are  wave 
structures.  From  the  upstream  region  to  the  location  slightly 
downstream  of  the  potential  core,  pressure  fluctuations  prop¬ 
agate  in  both  downstream  and  upstream  directions,  as  shown 
by  the  peak  correlation  lines  that  point  towards  both  direc¬ 
tions  as  the  delay  time  decreases.  However,  the  peak  correla¬ 
tion  lines  only  point  to  downstream  in  the  mixing  region, 
indicating  that  pressure  fluctuations  having  large  correlation 
values  propagate  downstream  in  this  region.  Although  the 
pressure  waves  at  the  reference  point  correlate  with  pressure 
fluctuations  along  the  lip-line  over  a  large  axial  distance,  the 
most  intense  region  of  correlation  lies  between  x/De  of  5  and 
10  or  roughly  from  the  third  shock  cell  to  the  end  of  the 
potential  core.  This  region  should  be  the  major  source  loca¬ 
tion  of  the  shock-associated  noise.  The  frequency  of  the  peak 
correlation  lines  pointing  towards  upstream  is  the  same  as 
that  of  the  screech  tone.  This  indicates  that  these  peak  corre¬ 


lation  regions  directly  reflect  the  correlation  of  the  screech 
tone  at  the  reference  point  with  the  pressure  fluctuations  on 
the  lip  line.  In  addition,  since  it  is  found  that  the  averaged 
propagation  velocity  of  these  streak  lines  is  around  O.Sa-^, 
the  pressure  waves  at  5De~  10De  on  the  lip  line  would  take 
an  amount  of  time  with  an  average  of  roughly  I  ()De/a^  to 
arrive  at  the  nozzle  exit.  It  can  be  seen  that  the  pressure 
waves  at  the  reference  point  are  indeed  strongly  correlated 
with  pressure  fluctuations  along  the  lip  line  from  5De-10De 
with  a  delay  time  centered  around  1  ()De/a^_.  Furthermore, 
the  pressure  waves  at  the  reference  point  are  strongly  corre¬ 
lated  with  the  pressure  fluctuations  near  the  nozzle  lip 
(x  =  0.0De)  at  a  later  time.  This  is  caused  by  the  helical  na¬ 
ture  of  the  screech  tone. 

III.  CONCLUSION  AND  SUMMARY 

Conical  C-D  nozzles  differ  from  traditional  CD  nozzles 
in  two  important  respects,  (1)  that  they  have  sharp  throats  and 
(2)  that  they  diverge  all  the  way  to  the  exit  creating  an  out¬ 
ward  component  of  velocity  in  the  exiting  flow.  These  two 
differences  each  cause  a  series  of  shock  diamonds  to  exist 
even  when  the  nozzle  is  operated  at  perfect  expansion,  when 
Mj  =  Md,  and  there  is  perfect  pressure  match  at  the  exit.  For 
all  conditions,  underexpanded  and  overexpanded  as  well  as 
perfectly  expanded,  the  presence  of  a  sharp  throat  creates  a 
second  set  of  shock  diamonds  independent  of  the  set  produced 
by  the  nozzle  exit  producing  a  double  diamond  pattern. 

The  lengths  of  shock  cells  from  conical  C-D  nozzles  are 
not  significantly  different  from  those  from  traditional  noz¬ 
zles.  The  shock  cell  length  is  an  important  parameter  in  con¬ 
trolling  the  frequency  of  both  broadband  shock-associated 
noise  and  screech,  so  conical  C-D  nozzles  produce  peak  fre¬ 
quencies  similar  to  those  of  traditional  C-D  nozzles.  The 
acoustic  character  distinctive  to  conical  C-D  nozzles  is  the 
presence  of  the  shock-associated  components  of  jet  noise 
even  when  the  nozzles  are  operated  at  perfect  expansion. 
Both  far-held  and  near-held  acoustic  measurements  on  coni¬ 
cal  C-D  nozzles  produce  the  same  shock-associated  noise 
signatures  one  would  expect  for  a  traditional  C-D  nozzle 
operated  off-design.  The  only  difference  being  that  these  sig¬ 
natures  exists  even  at  the  design  condition. 

Agreement  of  the  LES  employed  in  this  study  with  the 
measured  quantities  has  been  excellent.  Suitably  validated, 
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LES  provides  far  greater  detail  than  the  measurements  do 
and  provides  simultaneous  and  correlated  results  between 
pressure  and  velocity  over  the  entire  simulated  domain  while 
measurements  of  sound  and  flow  fields  are  taken  separately 
and  are  thus  uncorrelated. 

LES  reveal  the  behavior  inside  the  nozzle  explaining  the 
origin  of  the  second  set  of  shock  diamonds  and  the  central 
low-speed  region  in  the  wake  of  the  internal  Mach  disk.  In¬ 
stantaneous  LES  provides  a  clear  picture  of  the  pressure 
pockets  which  mark  the  large-scale  structures  as  they  con- 
vect  along  the  shear  layer.  The  temporal  pressure  fluctuations 
caused  by  the  passing  of  these  pressure  pockets  are  shown  to 
be  correlated  with  acoustic  pressure  fluctuations  at  a  refer¬ 
ence  location  outside  the  jet  in  the  region  through  which 
shock  associated  noise  is  known  to  pass.  The  high  correlation 
between  the  pressure  fluctuations  in  the  jet  shear  layer 
between  5  and  10  diameters  with  those  at  the  reference  loca¬ 
tion  are  convincing  evidence  that  this  is  the  source  location 
for  the  shock-associated  noise. 
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